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Motivated by the prospect of testing inflation from precision cosmic microwave background ob- 
servations, we present analytic results for scalar and tensor perturbations in single-field inflation 
models based on the application of uniform approximations. This technique is systematically im- 
provable, possesses controlled error bounds, and does not rely on assuming the slow-roll parameters 
to be constant. We provide closed-form expressions for the power spectra and the corresponding 
scalar and tensor spectral indices. 



Recent observations of the cosmic microwave back- 
ground (CMB) and of the distribution of galaxies have 
been largely consistent with inflation, not only with the 
result that the Universe is at critical density but also with 
the idea that the structure in the Universe arose from the 
gravitational collapse of adiabatic, Gaussian, and nearly- 
scale invariant primordial fluctuations, a key prediction 
of simple models of inflation [Q. The prospect of us- 
ing precision CMB observations to distinguish between 
models of inflation has focused attention on the predic- 
tive control of the theory underlying the generation of 
inflationary fluctuations. 

The form of the power spectrum alluded to above is 
predicted by models where inflation is caused by the dy- 
namics of a scalar field in a 'featureless' potential, evolv- 
ing in a friction-dominated 'slow-roll' regime. In more 
complex models (multiple fields, nontrivial potentials) it 
is possible to engineer spectra in a variety of ways; how- 
ever, present observations do not demand the considera- 
tion of dynamically exotic inflation models. 

Inflation predicts a spectrum of metric perturbations in 
the scalar (density) and tensor (gravitational wave) sec- 
tors, the vector component being naturally suppressed. 
Both scalar and tensor perturbations cause anisotropy in 
the CMB temperature Q . Tensor modes cause a poten- 
tially observable polarization of the CMB, a prime target 
for next-generation CMB observations || . 

In tests of the inflationary paradigm, a key point is 
that, in inflation, the scalar and tensor fluctuations are 
not independent. It is important to accurately quan- 
tify this relationship, e.g., via 'consistency relations' j|. 
However, approximations based on 'slow-roll' assump- 
tions lack error control and are not systematically im- 
provable ||. Attempts to cure this problem have so fai- 
led to restrictive assumptions and/or tedious mathemat- 
ical formulations. 

Our purpose in this Letter is to present an analysis of 
the inflationary perturbation spectrum by applying a uni- 
form approximation to the relevant mode equations. Our 
analysis does not make restrictive assumptions, possesses 
error bounds, and is systematically improvable. We also 
provide compact expressions for the scalar and tensor 
spectral indices n$ and tlt- Details will be reported else- 
where ||. 



The dynamical equations for gauge-invariant scalar 
perturbations can be found by linearizing the Einstein 
equations in the longitudinal gauge (see, e.g., Ref. Q). In 
terms of the standard choice of gauge invariant variable u 
(following the conventions of Ref. g), one of these equa- 
tions reduces to the compact form u" — Au — (z" / z)u = 
where z = a(j>' /h, h = a' /a, a is the scale factor and the 
prime denotes a derivative with respect to the conformal 
time r\ = J dt/a. This equation represents the dynamics 
of a free field with a time-dependent mass. The canonical 
form makes u the appropriate variable to quantize. 

To proceed further, we introduce the quantity C = u/ z, 
the intrinsic curvature perturbation of the comoving hy- 
persurfaces H . The dynamical equation for u forces ( to 
be nearly constant in the long wavelength limit k — > 0. 
This is true during the inflationary phase as well as in the 
post-reheating era. A computation of £ provides all the 
information needed (aside from the transfer functions) to 
extract the temperature anisotropy of the CMB. Details 
of this procedure can be found in Refs. |, @, 0. 

The calculation of power spectra involves computing 
two-point functions of quantum operators, e.g., 

/* (JjJ^ sin /ct" 

(0\u(r,, xjfifa, x + r)|0) = / -77-7— Pu(v,k). (1) 

The associated complex amplitude Uk(rf) satisfies 

u'i +(k 2 -^ju k = . (2) 

Solving Eqn. (^) is the fundamental problem in deter- 
mining the primordial power spectrum P u (or the power 
spectrum for £, here defined as Ps)- 

The case of the tensor perturbations proceeds analo- 
gously. The relevant mode equations are 

v'l + (V - ^ V = , (3) 

which correspond to the wave equation for a massless, 
minimally coupled, scalar field in an FRW Universe. 

One approach to analytically approximating Eqns. (|^) 
and (||) relies on the fact that exact solutions exist in the 
limits k 2 ^> \z"/z\ (short wavelength) and k 2 <C \z" / z\ 
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(long wavelength) or, as will be made more explicit below, 
as — krj — > oo and krj — > _ . For scalar perturbations, 

p -ikr) 

(k 2 ^\z"/z\, -fctj-Kjo), (4) 



u k 



A k z (k 2 « \z"/z\ , krj -> 0") . (5) 

Here, the short wavelength solution corresponds to the 
choice of an adiabatic vacuum for modes on length scales 
much smaller than the scale set by the curvature. The 
long wavelength solutions correspond to the growing 
mode on scales much larger than the Hubble length. Sim- 
ilar solutions exist for tensor perturbations. The problem 
of determining the power spectrum boils down to match- 
ing the asymptotic solutions, i.e., determining A k . 

To facilitate further analysis, it is useful to introduce 
the substitution 



Z "/Z E 

in Eqn. (Q), yielding 



(4(v) -i/4) h 2 



(6) 



u k ( V ) = 0. (7) 



In the standard slow-roll method, vs(j]) is taken to be 
constant in the intermediate regime — krj ~ 1 or k ~ aH 
(horizon-crossing for the fc-mode of interest), in which 
case the mode equations can be solved in terms of Bcsscl 
functions and these intermediate solutions matched to 
the short and long wavelength solutions || (see the dis- 
cussion in Ref. 0). Systematic improvement of the ap- 
proach runs into serious difficulties fj): as soon as vs{rj) 
is varying, the Bessel solution no longer holds. (A per- 
turbative improvement based on the leading order Bessel 
solution is given in Ref. See also, Ref. 

To overcome these difficulties, we will here apply the 
method of uniform approximation |]l3| . This method 
provides a single approximating solution which con- 
verges uniformly, for time-dependent v§, taking on the 
Liouville- Green (LG) or WKB form in both the short 
and long wavelength limits. The uniform approximation 
employs no intermediate matching, and, as opposed to 
the standard slow-roll approach, the approximation pro- 
cedure can be systematically improved and possesses an 
error control function || . 

We now turn to the specific form of the approximate 
solutions. Following Ref. jl3| we find a solution for u k {ifj 
containing a part valid to the left of the turning point, 
V = Vs (defined by k 2 = v s /fj 2 s ), and a part valid to the 
right. The leading order unnormalized solutions to the 
left (j] < fjs) and right {rj > fjs) of the turning point 
can be written in terms of Airy functions (Ai^ = Ai, 
Ai( 2 ) = Bi) as 



<< W 

> 



f<{v) 



9s(v) 



= /< /4 W5s 1/4 WAi (1 ' 2) [/<W], 

> > 



(8) 



^±2 J d ^^9s(v')] 1/2 j , (9) 



"si*)) 



(10) 



The general solution is a linear combination u k — Au^ + 
Bu k with the constants A and B ensuring that Ukiji) = 
(1/ '^/2k)e~ lkri in the limit k — » oo. In this limit, for 
well-behaved vs, f<(v) i s large and negative and we can 
employ the asymptotic forms 



{AiW,Ai^}(-x) = 



{cos, sin} (2x 3 / 2 /3- n/4) 



ri/a-i/4 



(11) 



Making the choices, A = \/n/2, B = we ob- 

tain the required adiabatic form of the solution at short 
wavelengths and as r\ — > — oo 14|. Note how the uni- 
form approximation makes the matching procedure very 
simple. 
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FIG. 1: Comparison of £ for a <f> chaotic inflation model 
calculated using the uniform approximation (crosses) and a 
numerical solution of the exact equations (solid line). 



The power spectra are computed in the limit 77 — > — . 
The I/77 2 pole dominates the behavior of the solutions 
and the Airy solution limits to the LG/WKB solution. 
The argument of the Airy functions, />(??), becomes 
large, allowing the use of the asymptotic forms 



{AiW.Ai^Kx) - 



-1/4 



exp 



{-,+} r V2 



{2,1}V^ 

(12) 

For computing the power spectra in the kr\ — > limit 
only the growing part of the solution is relevant: 



kr)—>0~ 



-iC 



n 



vs{r]) 



exp 



(13) 

where C is an irrelevant constant phase factor. The 
calculation for the tensor case proceeds in a similar 
fashion, with the replacement of us by vt defined by 
a" j a = [^xiv) ~ l/l] /v 2 - A numerical illustration of 
the accuracy of the approximation is provided in Fig. [I]. 

Once the approximate solutions to Eqns. (||) and 
(Q) have been found in the manner described above, 
the relevant power spectra can be easily computed via 
P s (k) = (k 3 /2TT 2 )\u k (r])/z\ 2 . Using Eqn.' (ph, we find 



Ps(k) 



1 



4tt 2 I z{ri) 1 2 v s {ri) 



exp 



vs 



(14) 
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The generalized spectral index for scalar perturbations 
can be obtained from the power spectrum via ns(k) = 
1 + d\nPs(k)/d\nk. It is important to note that the 
turning point is a function of k, since k = —vs{vs) Ms 
where vs{vs) is the value of vsiv) a t t ne turning point 
rj = fjs. Using this relation, one finds 



n s (k) = 4 + 2 



vs{rjs)k p drj' 



\Vs\ 



(15) 



The above analysis can be carried out for tensor per- 
turbations in an identical fashion with the replacement 
vs — > vt- The spectral index for tensors is given by 



n T (k) = 3 + 2 



v T (r] T )k 



drj 



y/grW) 



(16) 



The expressions (|l4|), (|T^), and ( ]T6| ) constitute our fun- 
damental result. 

In order to achieve a completely local expression, we 
can make additional approximations in the integrals in 
Eqns. and (|l(|). The integrand has a square-root 

singularity at the turning point while at the upper limit 
r\ goes to zero and the integrand vanishes linearly. There- 
fore, assuming vs{v) and vt(v) are well-behaved, we ex- 
pect the dominant contribution to arise from the lower 
limit. Assuming that vg and vt vary in such a way that 
one can expand them around the turning point in a Tay- 
lor series, to first order in derivatives, we find 

"Uv) - v 2 s {f)s) + 2v s {fjs)v' s {vs) (V ~ Vs) • (17) 



The integrals in Eqns. (|15|) and (|16|) can now be solved 
exactly, and the spectral indices are given by 



n s (k) 
n T (k) 



4 - 2v s (vs) 
3 — 2vt{t]t) 









1 - 




§) 














1 - 


('- 











v'sivs) , 

v 't(Vt] 
v T {r\T, 



Vt 



,(18) 
(19) 



As a simple test of the uniform approximation, we con- 
sider the case of constant vs and vt where exact results 
are available. In this case (|l4|), (|l5|), ( |l6| ) can be evalu- 
ated immediately. One finds 



Ps(k) 



p|*(fc) 1 

1 &vs 



H 

(10 
1 



724 



(-kr))- 2us+1 k 2 

(20) 
(21) 



ns = 4 — 2^5, riT = 3 — 2vt, 



where we have used Stirling's formula for the Gamma 
function in the exact power spectrum -P§ x (fc) (See, e.g., 
Ref. The spectral indices are found exactly, un- 

like in slow- roll expansions. All analytic examples so far 
used to test slow-roll results H p| , e.g., power-law in- 
flation, inflation near a maximum, and natural inflation 
have constant vg\ the uniform approximation recovers 
the spectral indices exactly in all of these cases. 



Error bounds follow from the computation of the vari- 
ation of the error control function 11131 



dv S 9 S 



-3/2 



9s 



:9s 



Vs) 2 



1 



t9s 



± 



24|/<|3/2' 



(22) 



This calculation allows the determination of the errors 
£1(77) and £2(77) in the two independent solutions and 

Uu (||). The error can be reduced if required by apply- 
ing higher-order uniform approximations. Details will be 
provided in Ref. ||. As an illustration, for the case of 
constant us, the absolute value of the relative error in the 
unnormalized solutions, over the full domain of interest, 
is bounded by 



£1,2 



< V2 



1 

6vs 



A 



724 



(23) 



where A ~ 1.04. The error bound has a weak subleading 
dependence on k arising from the end-points, falling as 
l/kf] far to the left of fjs and vanishing as k 2 rj 2 , in the 
limit 77 0~. The fact that the error bound is essentially 
/c-independent is consistent with the spectral indices be- 
ing exact in this case. The error in the power spectrum 
given in Eqn. ( pCf ) falls comfortably within the bound. 

It is instructive to recast our results in the more famil- 
iar language of the slow-roll expansion. The exact forms 
of z"/z and a" /a are 

2a 2 H 2 fl + e-^-^+^ + ^W) 



2a 2 H 2 I 1 - -e 



(25) 



where e = -H/H 2 , 5 = -$/(H(j>), and £2 = (e- 5)/H. 
The overdot represents a derivative with respect to the 
cosmic time t and H = a/ a. 

In the slow-roll approximation e, 5, and £2 are treated 
as small parameters, and one aims to solve Eqns. (j^) 
and (||) in terms of expansions in these parameters. The 
exact equations of motion for e and S are 

±=2e(e-S), A =2e(e (26) 

The conformal time rj can be written in terms of the slow- 
roll parameters via repeated integration by parts |l5|] : 



1 



1 



+ 2e(e-5) + 



(27) 



If e, S <C 1, £2 ~ 0(e 2 , S 2 ,eS), then, e and S are constant 
to leading order (and to this order, £ 2 — 0; the notation 
£2 makes explicit that £2 is a second-order quantity). At 
leading order, rj_^ -(1 + e)/aH [Cf. Eqn. @]. 

Our results (18) and ( |l9| ) correspond to a resummation 
of the slow-roll expansion to all orders. The procedure 
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of partially undoing this resummation to produce a less- 
accurate, but more familiar, expansion is as follows. We 
begin by substituting Eqn. ( pi| ) in the exact expression 
for Eqn. (|J). Expanding Eqn. (|27| ) to quadratic order 
in slow-roll parameters, a 2 H 2 ~ (1 + 2e + 7e 2 — 4e<5)/?7 2 , 
then yields, 



r,v' s ~ 2e(<5-e)-2£ 2 . 



(28) 



These second-order slow-roll results are new and can be 
compared against the results of Ref . jll) : both agree at 
leading order as well as for the case of power-law inflation, 
however the next-to-leading expressions differ in the gen- 
eral case. The second order approximations given above 
are consistent with the first order derivative expansions 
in Eqns. ( p"8| ) and ([l9|). The consistency of higher order 
derivative approximations can be determined by compar- 
ing against the computed error bounds at a given order 
of the uniform approximation. 



The desired form of the final result can be obtained by 
substituting these expressions in Eqn. (fill) 



n s {k) ~ l-4e + 25 + 2 



— - 7T ) e d 



, , 20 

-2 7T 

3 



3 2 



6)2, (29) 



where the subscripts on the slow-roll parameters de- 
note that they are evaluated at the turning point k = 
— vs(vs)/vs- Clearly this expression is of lower accuracy 
than (p"8"|). For example, one no longer recovers the exact 
results when vs is constant. The analogous result for the 
tensor spectral index is 

n T (k) ~ -2e - 2 (™ - n) e 2 Q + 2 f H - tt) e S . (30) 



To summarize, we have presented a controlled approx- 
imation for the computation of the inflationary pertur- 
bation spectrum which overcomes the matching problem 
of the conventional slow-roll approximation and system- 
atically improves the calculation in the generic case of 
varying v$ and vr. Future applications of this method 
include precision studies of generic predictions of inflation 
models JII|, reconstruction 15 , and studies of backreac- 
tion of quantum fluctuations on the background space- 
time metric. 
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